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Abstract. New over-ocean aerosol models are developed 
by integrating the inversion data from the Aerosol Robotic 
Network (AERONET) sun/sky radiometers with a database 
for the optical properties of tri-axial ellipsoid particles. The 
new aerosol models allow more accurate retrieval of aerosol 
optical depth (AOD) from the Moderate Resolution Imag- 
ing Spectroradiometer (MODIS) in the case of high AOD 
(AOD > 0.3). The aerosol models are categorized by us- 
ing the fine-mode fraction (FMF) at 550 nm and the single- 
scattering albedo (SSA) at 440 nm from the AERONET in- 
version data to include a variety of aerosol types found 
around the globe. For each aerosol model, the changes in the 
aerosol optical properties (AOPs) as functions of AOD are 
considered to better represent aerosol characteristics. Com- 
parisons of AODs between AERONET and MODIS for the 
period from 2003 to 2010 show that the use of the new 
aerosol models enhances the AOD accuracy with a Pearson 
coefficient of 0.93 and a regression slope of 0.99 compared 
to 0.92 and 0.85 calculated using the MODIS Collection 5 
data. Moreover, the percentage of data within an expected er- 
ror of ± (0.03 + 0.05 x AOD) is increased from 62 % to 64 % 
for overall data and from 39% to 51 % for AOD > 0.3. Er- 
rors in the retrieved AOD are further characterized with re- 
spect to the Angstrom exponent (AE), scattering angle (0), 


SSA, and air mass factor (AMF). Due to more realistic AOPs 
assumptions, the new algorithm generally reduces system- 
atic errors in the retrieved AODs compared with the current 
operational algorithm. In particular, the underestimation of 
fine-dominated AOD and the scattering angle dependence of 
dust-dominated AOD are significantly mitigated as results of 
the new algorithm’s improved treatment of aerosol size dis- 
tribution and dust particle nonsphericity. 


1 Introduction 

Aerosols exert a significant impact on climate change and air 

quality. The small airborne particles regulate the radiation 
budget through both direct and indirect effects (IPCC, 2007), 

specifically, by scattering and absorbing radiation and by 
modifying cloud microphysics. Aerosols are known to affect 
human health by causing and worsening respiratory illnesses 

(Pope and Dockery, 2006). Because the spatio-temporal dis- 
tribution of aerosols is highly variable, satellite observa- 
tions have been extensively utilized to quantify aerosol op- 
tical properties (AOPs) over wide areas and with fine spatio- 
temporal resolution. 
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Traditional 5-channel meteorological imagers, including 
single visible -band instruments aboard geostationary satel- 
lites, are used to continuously monitor aerosol optical depth 
(AOD), but have a limited ability to retrieve other parameters 
(Knapp et al., 2002; Wang et al., 2003; Kim et ah, 2008). In 
contrast, multi-spectral instruments on board low Earth orbit 
(LEO) satellites, such as the Advanced Very High Resolu- 
tion Radiometer (AVHRR), Sea-viewing Wide Field-of-view 
Sensor (SeaWiFS), and Moderate Resolution Imaging Spec- 
troradiometer (MODIS), can retrieve aerosol size informa- 
tion and absorptivity (Higurashi and Nakajima, 1999, 2002; 
Mishchenko et ah, 1999; Hsu et ah, 2004, 2006; Remer et ah, 
2005; Kim et ah, 2007; Levy et ah, 2007b). The Geostation- 
ary Ocean Color Imager (GOCI), which observes spectral ra- 
diances centered at 412, 443, 490, 555, 660, 680, 745, and 
865 nm from a geostationary orbit, has been used for hourly 
monitoring of AOD and to retrieve the fine-mode fraction 
(FMF) and aerosol types over East Asia (Lee et ah, 2010b). 

With their wide spatial and spectral coverage, the obser- 
vations made by the MODIS instruments aboard the Terra 
and Aqua satellites provide an unprecedented opportunity 
to infer AOPs. MODIS has 36 spectral bands ranging from 
0.41 to 15 pm with three different spatial resolutions (250 m, 
500 m, 1 km) and with 2300km-wide swath coverage. Since 
the launches of MODIS in 1999 for Terra and 2002 for Aqua, 
numerous efforts have been made to retrieve, evaluate, and 
improve the aerosol products obtained. The original opera- 
tional algorithms for dark vegetated areas (Kaufman et ah, 

1997) and oceans (Tanre et ah, 1997) are two very impor- 
tant additions. The MODIS algorithms have been frequently 
updated to improve the quality of retrieved data by modify- 
ing cloud-masking processes, aerosol models, and the sur- 
face reflectance database (Remer et ah, 2005; Levy et ah, 
2007a, b). Hsu et ah (2006) developed Deep-Blue algorithm, 
which is applicable to bright land surfaces including desert 
areas, to facilitate monitoring of dust aerosols over source 
regions. Consequently, the current MODIS operational algo- 
rithms provide the columnar aerosol amount (AOD) and size 
information (FMF, Angstrom exponent [AE]) for full cover- 
age of the Earth except for cloud- and snow-covered areas. 
In addition, the Deep-Blue algorithm also has the ability to 
retrieve the single-scattering albedo (SSA) of dust aerosols. 

MODIS aerosol products have been validated extensively 
to evaluate data quality. Preliminary comparisons of AOD 
from Terra-MODIS with that observed from Aerosol Robotic 
Network (AERONET) Sun/sky radiometers (Holben et ah, 

1998) showed that the AOD at 660 mn over the ocean dif- 
fered by only 2% on average, with negligible offset (Re- 
mer et ah, 2002), while the Terra-MODIS AOD over land 
was underestimated by about 14% except for coastal areas 
(Chu et ah, 2002). However, the comparison results for land 
varied significantly with location, partially due to different 
surface conditions and aerosol sources. On the contrary, a 
validation by Remer et ah (2008) showed an almost perfect 
regression slope for the AOD at 550 nm over land as calcu- 
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lated by AERONET and MODIS, but an underestimation of 
AOD over the ocean from Aqua-MODIS, in particular for 
high AODs. It should be noted that the land algorithm has 
been modified substantially to resolve better aerosol models 
and surface reflectance (Levy et ah, 2007a, b), whereas no 
substantial update has been made to the original ocean algo- 
rithm (Remer et ah, 2005, 2006). Therefore, the significant 
improvement in the regression slope over land is likely at- 
tributed to improved aerosol models to some extent, because 
the slope largely depends on data in the high AOD regime 
where the aerosol signal dominates the surface signal. 

In this paper, new aerosol models are introduced by inte- 
grating AERONET inversion data (Dubovik and King, 2000; 
Dubovik et ah, 2006) with single-scattering property data 
from a tri-axial ellipsoid database (Meng et ah, 2010). The 
AOPs of each aerosol model are used to calculate a lookup 
table (LUT) for spectral reflectances from MODIS. By us- 
ing the LUT, retrieval and validation of AODs are performed 
over the global ocean using Aqua-MODIS data to improve 
the underestimation of AOD reported by Remer et ah (2008). 
The validation results are compared with those from the cur- 
rent operational algorithm to further characterize the effects 
of the new aerosol models. 


2 MODIS ocean algorithms 

Remer et ah (2005, 2006) described the current MODIS Col- 
lection 5 operational over-ocean algorithm (hereafter, C005 
algorithm) in detail. The algorithm retrieves spectral AOD 
and FMF using spectral reflectances centered at 555, 650, 
860, 1240, 1630, and 2120 nm by comparing observed and 
pre -calculated reflectances. To this end, sophisticated for- 
ward radiative transfer simulations of the reflectances (i.e. 
LUT) need to be performed for various aerosol models, 
surface reflectances, and sun/satellite geometries. Because 
the top-of-atmosphere (TOA) reflectance consists of signals 
from both the surface and atmosphere, the algorithm spec- 
ifies surface reflectance in terms of the Fresnel reflection 
accounting for sea-surface roughness with a wind speed of 
6ms -1 and zero water-leaving radiance except at 550 nm 
where a water-leaving radiance of 0.005 is assumed. The at- 
mospheric contribution, aerosols in particular, is calculated 
using four fine-mode and five coarse-mode aerosols. Then, 
both spectral AOD and FMF (550 nm), a ratio of fine-mode 
AOD to total AOD, are simultaneously retrieved by mini- 
mizing the error between the observed and calculated re- 
flectances for each of the 20 combinations from the fine- and 
coarse-mode aerosol models with the fixed AOD retrieved 
from 860 nm. 

A schematic flowchart of the C005 algorithm and a test 
algorithm to evaluate the effects of the new aerosol models 
introduced in this paper is provided in Fig. 1. The test al- 
gorithm is designed to use the same spectral reflectances as 
the C005 algorithm in order to constrain other effects that 
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Fig. 1. Schematic flowchart of aerosol retrieval by the MODIS 
C005 algorithm (left column) and the test algorithm (right column). 
The test algorithm was designed to use the same observation data 
(“Mean_Reflectance_Ocean” in “MYD04” files) as the C005 algo- 
rithm to evaluate the effects of the new aerosol models only. The 
major difference between the two algorithms is the aerosol models 
used to calculate the LUT. 


can arise from different pre-processing of the data. We used 
the “Mean_Reflectance_Ocean” product in the “MYD04” 
files, which provides cloud- and sediment-masked mean re- 
flectance in 20 x 20 pixels of 500 m pixel-resolution data 
at seven wavelengths centered at 470, 555, 650, 860, 1240, 
1630, and 2120 nm. The product is the same as that used in 
the C005 algorithm. The major difference between the two 
algorithms is the aerosol model, and minor change is made 
to the inversion procedure. The test algorithm first retrieves 
AOD at 550 nm using every wavelength and aerosol model, 
and then selects the aerosol model that minimizes the stan- 
dard deviation of the seven different AODs retrieved from 
each wavelength. The final AOD is chosen according to the 
selected aerosol model. By doing so, each wavelength can 
contribute equally to selecting the aerosol model, whereas 
860 nm band has strong weighting in the C005 algorithm due 
to the calculation of spectral fitting error between observed- 
and simulated-reflectances with a perfect match at 860 nm in 
selecting an aerosol model. In addition, for the low AOD case 
(AOD < 0. 15), the test algorithm uses marine aerosol models 
and longer wavelengths (650, 860, 1240, 1630, 2120 mn) at 
which the water-body absorption is strong, whereas the al- 
gorithm uses full suite of aerosol models to be described in 
Sect. 3 and all the seven bands for AOD > 0.15. This proce- 
dure is done by first retrieving AOD using the marine aerosol 
models and the longer wavelengths over the whole areas and 
then retrieving again using the full suite of aerosol mod- 
els and all the wavelengths over the high AOD area deter- 
mined by the former procedure. The purpose of using differ- 


ent wavelength sets for the different aerosol loadings is to re- 
duce surface contribution to selecting an aerosol model par- 
ticularly for low AOD condition, and to use blue-wavelength 
which contains relatively high surface signal but also sensi- 
tivity to aerosol absorption for high AOD condition. There- 
fore, reflectance criterion instead of AOD can also be used 
to discriminate aerosol-signal-abundant condition. However, 
marine aerosol assumption for AOD <0.15 ties us to use 
AOD criterion even though there could be errors in discrim- 
inating high aerosol loading by the retrieved results. As a 
result, the algorithm simultaneously retrieves AOD and an 
aerosol model, i.e. FMF and SSA since the aerosol models 
are categorized by using these two parameters. Note that no 
combination between aerosol models is assumed in the test 
algorithm in contrast to the MODIS operational algorithm. 


3 New aerosol models 

Use of a radiative transfer model (RTM) to simulate satellite- 
observed TOA reflectance requires aerosol characteristics 
such as spectral refractive indices, size distribution, and non- 
sphericity to describe nonspherical particles. Otherwise, the 
spectral AOD, SSA, and phase function, which are derived 
from the aforementioned aerosol properties, are required. 
Thus, long-term AERONET inversion data that provides the 
AOPs for the globe can be used to simulate the satellite 
signal for various aerosol types. It should be noted that 
the AERONET observes ambient-columnar properties sim- 
ilar to those obtained from satellite observations, while in- 
situ measurements provide near-ground properties. More- 
over, the AERONET-retrieved AOPs represent the radiation 
field well for a wide scattering angle range, because the inver- 
sion data are retrieved to match the calculated radiation field 
with the observed sky radiances from the combined princi- 
pal/almucantar planes (Dubovik and King, 2000; Dubovik et 
al., 2006). The inversion data provide AOPs at 440, 675, 870, 
and 1020 nm, but for MODIS observations, the tri-axial ellip- 
soid database, introduced in the work of Meng et al. (2010), 
is required to expand the wavelength range up to 2120 nm. 

3.1 AERONET inversion data 

The quality-assured, “Level 2 Inversion All Points” data are 
used to derive aerosol models over the ocean for the test al- 
gorithm. To this end, data needs to be collected from specific 
AERONET stations chosen by distance from the ocean. The 
distance from the ocean is calculated by using geo-location 
information for each AERONET site and a high-resolution 
digital elevation model (DEM). The criterion for selecting 
the coastal stations was the distance within 7 km from the 
ocean. Figure 2 shows the 81 selected stations and the num- 
ber of inversion data available to date. Although large por- 
tions of the data are from the US and Europe, where an- 
thropogenic aerosols are dominant, the AERONET data in 
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Fig. 2. Global distribution of AERONET sun/sky radiometers lo- 
cated in coastal areas (81 stations) used to archive aerosol optical 
properties for the test algorithm. The colors represent the number of 
inversion data points at each site. AERONET stations within 7 km 
from the ocean were chosen as coastal stations. 


the downwind of North Africa are expected to provide AOPs 
of transported dust and biomass-burning aerosols. The data 
also cover dust aerosols transported from the Arabian Desert, 
both anthropogenic and dust aerosols over East Asia, and ma- 
rine aerosols over the remote ocean. However, only one site 
is located downwind of Southern Africa and a lack of data 
may lead to uncertainties in the AOPs of biomass-burning 
aerosols from the area. 

Aerosol models incorporated into satellite algorithms 
should account for the various aerosol types that exist and 
cause differences in the radiation field in order for the appro- 
priate aerosol models to be selected from the observed radi- 
ation field. Therefore, use of size and absorptivity to classify 
aerosol types is the most reasonable method for remote sens- 
ing applications, because the two parameters directly affect 
the radiation field (Dubovik et al., 2002; Levy et al., 2007a; 
Mielonen et al., 2009; Lee et al., 2010a). While MOD1S op- 
erational algorithms adopt fine- and coarse-mode aerosols 
separately and combine their signals during retrieval, the test 
algorithm adopts independent mixture-type models by clas- 
sifying aerosol types from the AERONET explicitly with re- 
spect to the FMF (550 nm) and the SSA (440 ran) (Lee et al., 
2010b). 

Figure 3 shows the number of aerosol events with spe- 
cific FMF and SSA values observed by AERONET through- 
out the globe and in coastal areas. The global data infers 
that aerosols from different locations have a wide range 
of FMF and SSA values, indicating the presence of vari- 
ous aerosol types with small to large particles and both ab- 
sorbing and non-absorbing. For FMF less than 0.4 (coarse- 
mode dominance), the SSA generally ranges from 0.85 
to 0.95, indicating absorption of blue-wavelengths by the 
coarse particle-dominated aerosols. Because Level 2 inver- 
sion data provides SSA for AOD (440 nm) > 0.4, the coarse 


particle-dominated aerosols mainly represent dust events. 
Note that sea salt, non-absorbing coarse-mode aerosols, 
generally occurs with low AOD values. For fine particle- 
dominated aerosols (FMF > 0.6), the range of values for SSA 
is wider than that for coarse particle-dominated aerosols. 
The high SSA values correspond to non-absorbing anthro- 
pogenic aerosols, such as sulfates and nitrates; whereas, the 
low SSA values imply the presence of black carbon (BC) 
(Hess et al., 1998; Wang and Martin, 2007). The major dif- 
ference between global and coastal data is highlighted by the 
lack of data in the extremely low SSA regime. The high rel- 
ative humidity (RH) in coastal areas, the aging of BC dur- 
ing transport, and the few AERONET stations in downwind 
of biomass-burning aerosols may cause the result. Both high 
RH and BC aging are known to increase SSA (Wang and 
Martin, 2007). 

Based on the representation of aerosol types classified us- 
ing FMF and SSA (Lee et al., 2010a), aerosol models are 
created by quantized square-bins over the FMF and SSA do- 
mains. Binning intervals of 0.1 and 0.05 are used for FMF 
and SSA, respectively. Each aerosol model is further catego- 
rized as a function of AOD by averaging AOPs between fore- 
and aft-medians of each AOD nodal point. If no data exist 
for a higher AOD bin, the AOPs of the previous bin is used. 
Consequently, the spectral AOD, SSA, and phase function 
are averaged over the three-dimensional domains of FMF, 
SSA, and AOD to be used as input data for the LUT cal- 
culations. Prior to the averaging, spectral AOD is normalized 
by itself at 550 nm and multiplied by each AOD nodal point. 
This method, however, has limited application to low AOD 
data, because SSA is retrieved only for AOD (440 nm) > 0.4. 
As an alternative, the SSA is assumed to be 0.99, regard- 
less of the wavelength, by considering the sea salt dominance 
in the low AOD regime over the ocean, while the other pa- 
rameters (spectral AOD, phase function) are compiled from 
AERONET inversion data. For marine aerosol models used 
for AOD < 0.15, only FMF criteria and AOD < 0.15 are used 
to calculate the AOPs, and SSA is assumed to be 0.99. As a 
result, a total of 23 aerosol models for AOD >0.15 and 9 
marine aerosol models for AOD <0.15 are created with the 
number of data points constrained to be greater than 10 for 
each aerosol model. The aerosol models cover FMF ranging 
from 0.2 to 1.0 for 0.85 < SSA < 0.95 (16 types), from 0.3 to 
1.0 for SSA > 0.95 (7 types), and from 0.1 to 1.0 (9 types) 
for marine aerosol models. Dimensions of the LUT are sum- 
marized in Table 1. Note that simple extrapolation is used 
to retrieve AOD of severe dust storms or smoke events with 
AOD higher than 3.6 (LUT limit). 

3.2 Tri-axial ellipsoid database 

The AERONET inversion data provide AOPs for wave- 
lengths ranging from 440 to 1020 nm, while the MOD1S 
observations cover the wavelengths from 470 to 2120 nm. 
To expand the wavelength range of AERONET AOPs, data 
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Fig. 3. The number of data points included in each FMF (550 nm) and SSA (440 nm) bin, archived from the AERONET inversion data over 
the globe (left) and coastal areas (right). The data were sorted into intervals of 0.05 and 0.01 for FMF and SSA, respectively. The AERONET 
stations in the coastal area are shown in Fig. 2. 


Table 1 . LUT dimensions for the MODIS over-ocean algorithm. 


Variable Name 

No. of 
Entries 

Entries 

Wavelength (A.) 

7 

470, 555, 650, 860, 1240, 
1630, 2120 nm 
(band 3, 4, 1,2, 5, 6, 7, 
respectively) 

SZA (6> 0 ) 

8 

0, 10, ..., 70° 

SAZA (0 S ) 

8 

0, 10, ..., 70° 

RAA ((/>) 

19 

0, 10, ..., 180° 

AOD (r) 

9 

0.0, 0.1, 0.3, 0.6, 1.0, 1.5, 
2.1, 2.8, 3.6 

Aerosol Model 

23 

Classified by FMF and SSA 
from AERONET inversion 
data 


SZA: solar zenith angle, SAZA: satellite zenith angle, RAA: relative azimuth 
angle. 


from the tri-axial ellipsoid database (Meng et al., 2010) 
are used in this study. The database containing the single- 
scattering properties of individual tri-axial ellipsoidal parti- 
cles was computed using the Lorentz-Mie code (Bohren and 
Huffman, 1983), the T-matrix code (Mishchenko and Travis, 
1998), the Amsterdam discrete dipole approximation (DDA) 
code (Yurkin and Hoekstra, 2009), and the improved geo- 
metric optics method (1GOM) code (Yang and Liou, 1996; 
Yang et al., 2007; Bi et al., 2009). Because tri-axial ellip- 
soidal shapes include spheres and ellipsoids with a number 
of aspect ratios, the application of the database can be ex- 
panded to non-dust aerosols. The database provides extinc- 
tion efficiency, SSA, phase matrix, etc., for various refractive 
indices, size parameters, and aspect ratios, thus the proxy of 


the AERONET AOPs (spectral AOD, SSA, and phase func- 
tion) can be extracted when given the spectral refractive in- 
dices, size distribution, aspect ratios, and nonsphericity. 

The procedure to calculate the AOPs of the predefined 
aerosol models at the longer wavelengths (1 > 1240 nm) is 
summarized as follows. 

1. Extracting extinction efficiency, SSA, and phase func- 
tion for individual spherical/nonspherical particles from 
the database by inputting refractive indices, size param- 
eter, and aspect ratios. Note that the RTM we used re- 
quires spectral AOD, SSA, and phase function to cal- 
culate TOA reflectances. The refractive indices for each 
aerosol model are from the MODIS operational algo- 
rithm since AERONET inversion data does not provide 
the values at the longer wavelengths. 

2. Integrating the AOPs of individual particles using pre- 
defined size distributions by assuming spherical par- 
ticles. The size distribution of each aerosol model is 
calculated by averaging AERONET inversion data the 
same as AOPs for the shorter wavelengths. 

3. Doing the same procedure as (2) for pure nonspheri- 
cal particles. Thus, integrating the AOPs of individual 
particles in both size distribution and aspect ratio dis- 
tribution space. The same aspect ratio distribution as 
AERONET inversion algorithm is used for consistency 
with the AOPs calculated for the shorter wavelengths. 

4. Calculating weighting mean of the AOPs for pure 
spheres and pure nonspherical particles using %spheric- 
ity. The mean values of %sphericity of each aerosol 
model are calculated by averaging AERONET inversion 
data. 
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Fig. 4. TOA reflectance difference at 860 nm between MODIS aerosol models and new aerosol models with respect to AOD for a given 
geometry. SZA (9 0 ) and SAZA (0 S ) are assumed to be 50° for both tests and RAA ( <p ) is assumed to be 90° (left) and 170° (right), resulting 
in scattering angle (0) of 114° and 172°, respectively. Corresponding MODIS aerosol models to the new aerosol models are created by 
combining F2 (“Water Soluble”) and C8 (“Dust-like type”) aerosol models from Remer et al. (2006) for SSA < 0.95 and F4 (“Water Solublt 
with humidity”) and C8 for SSA > 0.95 using FMF values from new aerosol models. Negative values are shown in line-fill. 


More specifically, the refractive indices of “water solu- 
ble” and “water soluble with humidity” from the C005 
aerosol models (Remer et al., 2006) are used for fine-mode 
(0.85 < SSA < 0.95 and SSA > 0.95, respectively), and the 
refractive indices of “dust-like type” and “wet sea salt type” 
are used for coarse-mode (AOD >0.1 and AOD = 0.1, re- 
spectively). Although the refractive indices of the MODIS 
aerosol models are not completely consistent with those of 
the new aerosol models, the low sensitivity of satellite signals 
to the fine-mode aerosols at longer wavelengths (> 1240 nm) 
due to their low AOD (high AE) and the relatively well- 
known absorption properties of dust (almost non-absorbing) 
(Hsu et al., 2004; Remer et al., 2005) are expected to result in 
smaller errors compared to excluding the longer wavelengths 
in the retrieval algorithm. A test retrieval using only shorter 
wavelengths (not shown for brevity) showed that exclusion 
of the longer wavelengths resulted in larger uncertainties than 
its counterpart particularly for determining aerosol size infor- 
mation. Nonetheless, more accurate optical property infor- 
mation for the longer wavelengths is still desirable for further 
improvement in the aerosol models and thus retrieval results. 

With regard to nonsphericity and aspect ratios, mean “% 
sphericity” from the AERONET data and fixed spheroid 
mixture distribution are used as described in Dubovik 
et al. (2006). The dataset compiled by combining the 
AERONET inversion data and tri-axial ellipsoid database 
data is used as input data into RTM calculations. The op- 
tical properties of the 23 aerosol models are summarized 
in Table 2. To calculate the LUT, a discrete ordinate radia- 
tive transfer (DISORT) code implemented into the libRadtran 
software package (Mayer and Kylling, 2005) is used. Due to 
the advantage of directly inputting the phase function into 
the libRadtran package, nonsphericity, which mainly affects 
phase function, can be readily handled with the software. 
However, the polarization effect is not included the same as 
in the C005 algorithm. 


4 Sensitivity study 

The possible differences in retrieved AOD between the C005 
algorithm and test algorithm are analysed before applying the 
developed algorithm to actual TOA reflectance data. Figure 4 
shows the difference in the calculated TOA reflectance values 
between the C005 aerosol models and the new aerosol mod- 
els for different AOD values at 860 nm, the reference wave- 
length used in selecting aerosol models from the C005 algo- 
rithm. The other wavelengths show similar tendencies. The 
corresponding C005 aerosol models to the new aerosol mod- 
els are created by combining the F2 (“water soluble”) and 
C8 (“dust-like type”) aerosol models in Remer et al. (2006) 
for SSA < 0.95 and the F4 (“water solublt with humidity”) 
and C8 for SSA > 0.95 by using FMF values from the new 
aerosol models. The results generally show that the C005 
aerosol models overestimate the TOA reflectance, i.e., under- 
estimation of AOD with a given TOA reflectance, compared 
to the new aerosol models, and the overestimation tends to 
increase with AOD and absorptivity. The increasing overes- 
timation with AOD is mainly due to the increasing particle 
size with AOD from the new aerosol models. Increasing par- 
ticle size results in decreasing TOA reflectance due to a de- 
crease in the back-scattering fraction of reflected radiation. 
Only slight underestimations are found in the dust-dominated 
models (H2, H3, H4, M2) and in the non-absorbing models 
(N1 through N7) in the side-scattering case for low AOD. 
On the other hand, consistent overestimation is present in 
the back-scattering direction, indicating difference in phase 
function between spherical and nonspherical particles. The 
phase function of spherical particle has lower values at side- 
scattering angles and higher values at back-scattering angles 
compared to their nonspherical counterparts. The underes- 
timation from the C005 aerosol model for non-absorbing 
fine-mode (N5 through N7) for low AOD is partially due to 
larger particle size of the F4 aerosol model compared with 
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Table 2. Spectral SSA (upper) and asymmetry factor (lower) with respect to wavelength, Angstrom exponent (AE), effective radius (r e ff), 
and sphericity of the new aerosol models used for the test algorithm. H models (0.85 < SSA < 0.90) and M models (0.90 < SSA < 0.95) 
cover FMF ranging from 0.2 to 1.0, while N models (SSA >0.95) cover from 0.3 to 1.0. The minimum and maximum values are shown 
because of AOD dependence. The optical properties are interpolated to spectral response functions of MODIS bands in RTM calculations. 


Model No. 
(SSA, FMF) 

440 nm 

675 nm 

870 nm 

1 020 nm 

1240 nm 

1 640 nm 

2120 nm 

AE 

reff (l*m) 

sphericity 

HI 

0.875-0.879 

0.932-0.933 

0.940-0.946 

0.944-0.951 

0.990-0.991 

0.985-0.986 

0.994-0.995 

0.16-0.22 

0.76-1.08 

i-i % 

(0.875, 0.25) 

0.752-0.763 

0.723-0.733 

0.728-0.739 

0.732-0.742 

0.729-0.731 

0.733-0.734 

0.736-0.738 




H2 

0.877-0.877 

0.920-0.922 

0.929-0.932 

0.934-0.938 

0.989-0.989 

0.982-0.983 

0.993-0.993 

0.27-0.28 

0.61-0.63 

2-3 % 

(0.875, 0.35) 

0.719-0.726 

0.696-0.698 

0.705-0.710 

0.716-0.720 

0.725-0.726 

0.731-0.732 

0.735-0.735 




H3 

0.874-0.884 

0.902-0.914 

0.916-0.925 

0.920-0.928 

0.988-0.988 

0.979-0.979 

0.992-0.992 

0.32-0.43 

0.53-0.54 

9-10% 

(0.875, 0.45) 

0.701-0.704 

0.674-0.679 

0.673-0.681 

0.690-0.699 

0.718-0.718 

0.727-0.728 

0.730-0.731 




H4 

0.877-0.880 

0.883-0.902 

0.898-0.913 

0.905-0.916 

0.986-0.986 

0.974-0.975 

0.991-0.991 

0.43-0.60 

0.47-0.50 

17-21% 

(0.875, 0.55) 

0.695-0.702 

0.664-0.667 

0.657-0.661 

0.674-0.675 

0.702-0.709 

0.715-0.720 

0.723-0.725 




H5 

0.873-0.873 

0.871-0.871 

0.881-0.881 

0.882-0.882 

0.984-0.984 

0.969-0.969 

0.989-0.989 

0.79-0.79 

0.45-0.45 

28-28 % 

(0.875, 0.65) 

0.702-0.702 

0.659-0.659 

0.651-0.651 

0.654-0.654 

0.684-0.684 

0.704-0.704 

0.717-0.717 




H6 

0.875-0.882 

0.860-0.873 

0.855-0.865 

0.854-0.857 

0.981-0.982 

0.951-0.955 

0.982-0.984 

1.24-1.28 

0.34-0.36 

67-74% 

(0.875, 0.75) 

0.706-0.736 

0.643-0.659 

0.627-0.63 1 

0.627-0.634 

0.617-0.639 

0.642-0.665 

0.675-0.690 




H7 

0.874-0.880 

0.846-0.876 

0.827-0.868 

0.814-0.861 

0.975-0.979 

0.933-0.944 

0.974-0.978 

1.33-1.48 

0.28-0.32 

74-74% 

(0.875, 0.85) 

0.703-0.735 

0.618-0.663 

0.584-0.617 

0.580-0.614 

0.591-0.602 

0.622-0.623 

0.663-0.666 




H8 

0.877-0.881 

0.850-0.869 

0.814-0.858 

0.787-0.842 

0.966-0.969 

0.891-0.899 

0.947-0.949 

1.65-1.69 

0.21-0.23 

96-97 % 

(0.875, 0.95) 

0.704-0.713 

0.601-0.630 

0.547-0.575 

0.523-0.546 

0.517-0.536 

0.526-0.530 

0.585-0.594 




Ml 

0.909-0.912 

0.956-0.963 

0.962-0.971 

0.965-0.973 

0.991-0.991 

0.985-0.985 

0.994-0.994 

0.18-0.19 

0.70-0.78 

2-3 % 

(0.925, 0.25) 

0.746-0.752 

0.723-0.730 

0.726-0.734 

0.732-0.741 

0.727-0.729 

0.732-0.733 

0.736-0.738 




M2 

0.912-0.914 

0.940-0.950 

0.950-0.959 

0.955-0.961 

0.989-0.990 

0.982-0.984 

0.993-0.994 

0.23-0.36 

0.62-0.76 

2-6% 

(0.925, 0.35) 

0.718-0.723 

0.690-0.703 

0.694-0.707 

0.706-0.715 

0.724-0.727 

0.731-0.734 

0.734-0.737 




M3 

0.916-0.916 

0.932-0.933 

0.941-0.942 

0.946-0.948 

0.987-0.988 

0.978-0.980 

0.992-0.992 

0.52-0.56 

0.52-0.60 

8-10% 

(0.925, 0.45) 

0.707-0.715 

0.681-0.687 

0.686-0.688 

0.698-0.700 

0.713-0.716 

0.725-0.727 

0.731-0.732 




M4 

0.922-0.923 

0.929-0.930 

0.935-0.936 

0.940-0.940 

0.985-0.985 

0.974-0.974 

0.990-0.990 

0.74-0.76 

0.45-0.47 

12-16% 

(0.925, 0.55) 

0.707-0.712 

0.670-0.673 

0.670-0.672 

0.683-0.685 

0.696-0.697 

0.714-0.715 

0.725-0.726 




M5 

0.922-0.922 

0.921-0.924 

0.925-0.927 

0.928-0.930 

0.983-0.983 

0.968-0.969 

0.988-0.988 

0.91-1.02 

0.38-0.43 

19-28% 

(0.925, 0.65) 

0.706-0.714 

0.656-0.660 

0.650-0.650 

0.659-0.661 

0.672-0.677 

0.696-0.701 

0.714-0.718 




M6 

0.923-0.924 

0.909-0.924 

0.905-0.924 

0.905-0.925 

0.980-0.980 

0.956-0.960 

0.984-0.985 

1.21-1.28 

0.34-0.35 

28—40% 

(0.925, 0.75) 

0.702-0.703 

0.639-0.641 

0.622-0.628 

0.629-0.633 

0.639-0.650 

0.671-0.678 

0.699-0.704 




M7 

0.929-0.929 

0.909-0.918 

0.896-0.909 

0.889-0.903 

0.974-0.977 

0.932-0.943 

0.973-0.977 

1.38-1.53 

0.28-0.31 

54-71 % 

(0.925, 0.85) 

0.717-0.723 

0.637-0.652 

0.601-0.619 

0.595-0.611 

0.579-0.603 

0.611-0.631 

0.662-0.673 




M8 

0.929-0.934 

0.909-0.937 

0.888-0.934 

0.871-0.930 

0.968-0.970 

0.899-0.908 

0.950-0.951 

1.48-1.70 

0.23-0.26 

69-92 % 

(0.925, 0.95) 

0.721-0.729 

0.639-0.660 

0.584-0.614 

0.557-0.589 

0.528-0.561 

0.526-0.540 

0.580-0.587 




N1 

0.949-0.959 

0.978-0.981 

0.980-0.983 

0.982-0.985 

0.991-0.992 

0.985-0.987 

0.994-0.994 

0.32-0.38 

0.62-0.71 

13-15% 

(0.975, 0.35) 

0.717-0.719 

0.698-0.698 

0.701-0.702 

0.709-0.709 

0.725-0.726 

0.730-0.731 

0.731-0.732 




N2 

0.948-0.959 

0.970-0.976 

0.976-0.981 

0.978-0.983 

0.990-0.990 

0.984-0.984 

0.993-0.993 

0.49-0.57 

0.58-0.62 

7-16% 

(0.975, 0.45) 

0.710-0.732 

0.686-0.698 

0.688-0.696 

0.697-0.705 

0.713-0.722 

0.722-0.730 

0.728-0.733 




N3 

0.953-0.961 

0.960-0.970 

0.964-0.972 

0.965-0.974 

0.988-0.988 

0.979-0.981 

0.990-0.991 

0.75-0.81 

0.42-0.45 

10-19% 

(0.975, 0.55) 

0.701-0.730 

0.664-0.684 

0.667-0.676 

0.680-0.685 

0.696-0.708 

0.708-0.720 

0.718-0.726 




N4 

0.959-0.963 

0.955-0.963 

0.952-0.965 

0.952-0.966 

0.986-0.987 

0.975-0.978 

0.988-0.989 

0.90-1.01 

0.36-0.41 

19-20% 

(0.975, 0.65) 

0.717-0.728 

0.660-0.674 

0.649-0.661 

0.658-0.666 

0.679-0.684 

0.699-0.702 

0.715-0.717 




N5 

0.961-0.963 

0.957-0.958 

0.953-0.955 

0.952-0.954 

0.983-0.984 

0.968-0.971 

0.985-0.985 

1.17-1.26 

0.32-0.35 

24-38 % 

(0.975, 0.75) 

0.717-0.732 

0.649-0.665 

0.629-0.644 

0.634-0.646 

0.648-0.656 

0.674-0.680 

0.699-0.706 




N6 

0.964-0.968 

0.959-0.961 

0.954-0.956 

0.951-0.953 

0.981-0.982 

0.957-0.958 

0.978-0.979 

1.39-1.49 

0.28-0.31 

48-68 % 

(0.975, 0.85) 

0.730-0.744 

0.652-0.671 

0.617-0.635 

0.610-0.624 

0.598-0.612 

0.625-0.630 

0.665-0.668 




N7 

0.969-0.970 

0.963-0.968 

0.956-0.963 

0.949-0.959 

0.978-0.980 

0.936-0.943 

0.961-0.963 

1.51-1.73 

0.24-0.28 

78-91 % 

(0.975, 0.95) 

0.734-0.755 

0.653-0.699 

0.600-0.654 

0.573-0.626 

0.542-0.584 

0.536-0.546 

0.565-0.593 





the fine-mode of the new aerosol model counterparts for low 
AOD. As a result, general underestimation of AOD is ex- 
pected from the C005 algorithm compared with the test al- 
gorithm, and this may have caused the underestimation of 
high AOD in the validation results represented in Renter et 
al. (2008). 

A redundancy test is performed to ensure whether all the 
aerosol models result in distinct spectral reflectance feature 
to be selected by the retrieval procedure. The test is carried 
out by retrieving AOPs, using calculated TOA reflectances 


as a proxy of observation data for the whole LUT dimen- 
sions shown in Table 1 (total 283 176 data points). Two dif- 
ferent proxy data are created by adding maximum random 
error of 3% and 10%. Figure 5 compares the exact solu- 
tions and the retrieved results. The data points with regard to 
various geometries and the other variables, except the tar- 
get variable, are averaged at each nodal point and shown 
with one-standard deviation intervals. The comparison re- 
sults show that each aerosol model (i.e. FMF and SSA) can 
be retrieved with small errors for the 3 % error case, while 
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Fig. 5. Comparison between input variables (AOD, FMF, SSA) and retrieved variables from the present algorithm in LUT space. Tests are 
performed with synthesized data including maximum random error of 3 % (upper) and 1 0 % (lower). Mean and standard deviation values are 
shown for each calculation point for LUT. 


significant misselection occurs for the 10% error case. The 
wrong retrieval is more prominent for SSA than FMF, show- 
ing lower accuracy in retrieving absorptivity than size in- 
formation from MOD1S. However, all the aerosol models 
are still required for the retrieval since highly-absorbing and 
non-absorbing aerosols are still discernable in some degree 
even for 10% error case. 

5 Results and evaluation 

The effects of the new aerosol models on AOD retrieval 
are evaluated by comparing the AODs from AERONET 
and MODIS data retrieved by the C005 and the test algo- 
rithm. Eight years of spectral reflectance data (2003-2010) 
observed from Aqua-MODIS are collected and processed 
to retrieve AOD using the new aerosol models. The AOD 
data from the C005 algorithm are also processed to compare 
with the AERONET observations, thereby allowing valida- 
tion results from both algorithms to be compared. In this 
study, the “Effective_Optical_Depth_Average_Ocean” data in 
the MYD04 files are used for the C005 algorithm. Overall 


statistical scores and systematic errors are compared to char- 
acterize various error sources. 

5.1 Overall evaluation 

Figure 6 compares AODs between AERONET and MODIS 
over the global ocean from 2003 to 2010. Three different re- 
sults retrieved from the C005 algorithm and the test algo- 
rithm using two different inversion procedures described in 
Sect. 2 are shown to investigate effects of new aerosol mod- 
els and inversion methods on AOD retrieval accuracy. For 
this comparison, collocation was made within ± 30 min in 
time and 25 km in space similar to the method proposed by 
Ichoku et al. (2002). However, we sampled satellite pixels 
by calculating actual distances between each AERONET lo- 
cation and the MODIS pixels for both test and operational 
datasets, while Ichoku et al. (2002) selected 5x5 MODIS 
pixels with an AERONET station located in the middle of 
the grid regardless of viewing angle. A criterion is applied 
for the number of data points, which requires at least 5 and 
2 data points for MODIS and AERONET, respectively. Be- 
cause the test algorithm tended to retrieve more data than the 
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Fig. 6. Comparison of AOD between AERONET and MODIS over the global ocean for the period from 2003 to 2010. The MODIS AODs are 
from the C005 algorithm (left) and the test algorithm (middle and right) with new aerosol models. Two different inversion procedures using 
standard deviation of spectral AOD (middle) and MODIS operational inversion (right) are applied for the test algorithm. The collocation 
criteria of ± 30 minutes in time and 25 km in space were used. The gray dots represent all data points, whereas black dots with one- 
standard deviation interval represent mean AODs in 20 equal-number-of-data bins with respect to the AERONET data. The solid line is 
from the regression equation, while the dotted and dashed lines are the one-to-one line and the MODIS expected error (EE) line showing 
± (0.03 + 0.05 x AOD), respectively. Only data points overlapping between the two algorithms are compared. Originally, the number of data 
points was 3106 for the C005 algorithm and 3578 for the test algorithm. The statistics shown are the Pearson coefficient (R), root mean 
squared error (RMSE), mean bias (MB), and the number of data points (IV). 


Table 3. Statistics for a comparison of AOD retrieved from the MODIS C005 algorithm and AERONET observations from 2003 to 2010 
over the global ocean. The numbers in the parentheses are for AOD (AERONET) > 0.3. 


MODIS 

OP 

2003 

2004 

2005 

2006 

2007 

2008 

2009 

2010 

Overall 

R 

0.91 

0.93 

0.92 

0.87 

0.94 

0.92 

0.92 

0.94 

0.92 

Slope 

0.78 

0.87 

0.82 

0.77 

0.9 

0.9 

0.88 

0.87 

0.85 

y-intercept 

0.04 

0.02 

0.03 

0.04 

0.02 

0.03 

0.03 

0.03 

0.03 

Percentage 

66% 

62% 

62% 

62% 

58% 

62% 

64% 

65% 

62% 

within EE 

(38%) 

(37%) 

(40%) 

(40%) 

(28%) 

(43 %) 

(44%) 

(54%) 

(39%) 

RMSE 

0.06 

0.07 

0.05 

0.05 

0.06 

0.06 

0.05 

0.05 

0.06 


(0.11) 

(0.13) 

(0.09) 

(0.10) 

(0.11) 

(0.10) 

(0.07) 

(0.07) 

(0.11) 

MB 

0.00 

0.00 

0.00 

0.00 

0.01 

0.01 

0.02 

0.01 

0.01 


(-0.06) 

(-0.05) 

(-0.06) 

(-0.06) 

(-0.02) 

(-0.01) 

(-0.02) 

(-0.04) 

(-0.04) 

N 

393 

451 

441 

346 

377 

336 

347 

264 

2955 


(50) 

(77) 

(47) 

(30) 

(38) 

(30) 

(29) 

(24) 

(325) 


C005 algorithm due to the absence of a quality-control proce- 
dure used in C005 algorithm, only overlapping data retrieved 
by both algorithms were used for quantitative comparisons. 

The validation results show that the AOD data from the 
C005 algorithm are highly correlated with the observations, 
but, on average, tend to be underestimated, with a Pearson 
coefficient of 0.92 and a regression slope of 0.85. The neg- 
ative bias of the slope is caused by overestimation in the 
low AOD regime (AOD < 0.2) and underestimation in the 
high AOD regime (AOD > 0.3). Meanwhile, the new aerosol 
models improve the slope significantly (0.99-1.02) with a 
comparable correlation coefficient (0.93) regardless of the 
inversion methods. Only small differences are observed be- 
tween the two results, showing a slight increase in statis- 


tics for the new inversion compared with the C005 inver- 
sion. However, the new inversion may not guarantee im- 
proved results for the C005 algorithm because of differences 
in the LUT. Since the difference is negligible between the 
two methods, only results from the new inversion will be 
shown for further analyses. From the statistics summarized 
in Tables 3 and 4, the aerosol models clearly improve almost 
all the statistics analysed in this study. The slope and the 
percentage of data within an expected error improved from 
62 % to 64 % overall and from 39 % to 51 % for AOD > 0.3. 
These improvements are particularly noticeable for the high 
AOD regime at which the aerosol signal dominates the other 
contributions, such as Rayleigh scattering and surface re- 
flectance. The major reason for the improvement in AOD 
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Table 4. Same as in Table 3 except for the test algorithm. 


This study 

2003 

2004 

2005 

2006 

2007 

2008 

2009 

2010 

Overall 

R 

0.94 

0.96 

0.93 

0.88 

0.94 

0.92 

0.93 

0.94 

0.93 

Slope 

1.01 

1.00 

0.97 

0.91 

1.01 

1.00 

1.01 

1.00 

0.99 

y-intercept 

0.01 

0.01 

0.01 

0.02 

0.01 

0.01 

0.01 

0.00 

0.01 

Percentage 

64% 

65% 

65% 

62% 

59% 

61% 

66% 

67% 

64% 

within EE 

(50%) 

(55%) 

(48%) 

(43 %) 

(36%) 

(60%) 

(51%) 

(70%) 

(51%) 

RMSE 

0.05 

0.05 

0.05 

0.06 

0.06 

0.06 

0.05 

0.05 

0.05 


(0.08) 

(0.08) 

(0.08) 

(0.10) 

(0.12) 

(0.10) 

(0.09) 

(0.07) 

(0.09) 

MB 

0.01 

0.01 

0.00 

0.00 

0.01 

0.01 

0.01 

0.00 

0.01 


(0.02) 

(0.01) 

(0.00) 

(-0.01) 

(0.05) 

(0.04) 

(0.02) 

(0.01) 

(0.02) 

N 

393 

451 

441 

346 

377 

336 

347 

264 
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Fig. 7. AE dependence of retrieval errors for the C005 algorithm (left) and the test algorithm (right) for overall data (upper) and data for 
AOD > 0.3 (lower). The data are sorted in 20 and 10 equal-number-of-data bins for the overall case and high AOD case, respectively. The 
dots and bars represent mean and one-standard deviation intervals of the retrieval errors, respectively, while the triangles represent the mean 
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retrieval is the consideration of absorbing fine-mode aerosols 
and changing the size distribution as a function of AOD. Both 
factors are expected to increase AOD for a given TOA re- 
flectance as shown in Fig. 4 and other previous studies (Levy 
et al., 2007b; Wang and Martin, 2007; Jethva et al., 2010). 
Note that the current MODIS algorithm adopts four different 
water-soluble aerosol models with fixed radii insensitive to 
AOD for the fine-mode cluster. 


5.2 Error characteristics 

Errors in AOD can arise from various sources including in- 
correct assumptions about surface reflectance and aerosol 
type, status of sensor calibration, observation geometry, 
etc. In addition, specific observation environments can bias 
results. Levy et al. (2010) evaluated the C005 over-land 
AOD data with regard to AE, cloud fraction, surface type 
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Fig. 8. Scattering angle dependence of retrieval errors for the C005 algorithm (left) and the test algorithm (right) for overall data (upper) and 
data for AOD > 0.3 (lower). 


characteristics, and observation geometry to characterize 
systematic error sources. Thus, the analyses to be presented 
here together with the work of Levy et al. (2010) represent a 
complete evaluation of the MODIS data retrieved over land 
and ocean. 

Figure 7 shows the AE dependence of the retrieval er- 
rors for the C005 and the test algorithms. The data were 
sorted into 20 and 10 equal-number-of-data bins for the over- 
all data and AOD > 0.3, respectively. As shown in the result, 
the C005 algorithm tends to overestimate coarse-dominated 
AOD (AE < 0.8) and underestimate fine-dominated AOD 
(AE>1.6). The underestimation of fine-dominated AOD 
worsens for high AOD cases, suggesting that the underes- 
timation of high AOD from the C005 algorithm shown in 
Fig. 6 is mainly caused by fine-dominated cases, while reli- 
able retrievals are perfonned for dust aerosols. The retrieved 
AODs are less stable (stability inferred by the magnitude 
of one standard deviation interval) for fine-dominated AOD 
than coarse-dominated AOD, while stable retrieval is ob- 
served for 0.75 < AE <1.4 partially due to the relatively low 
AOD. In the case of AOD > 0.3, although there are system- 
atic underestimations, the algorithm shows small mean bi- 
ases (MBs) for strong dust events (AE < 0.4) with a higher 
stability than the fine-dominated case (AE > 1.4). 


The AE dependence of MB for the test algorithm is re- 
duced overall compared with the C005 algorithm, but the 
test algorithm still has a tendency to overestimate coarse- 
dominated AOD and to underestimate fine-dominated AOD. 
The AE dependence is reduced further for AOD >0.3, 
but the AODs are distinctly overestimated for a coarse- 
dominated case (AE < 0.3). The improvement in the fine- 
dominated AOD for a high AOD case can be explained by 
the inclusion of the absorbing aerosol models and size dis- 
tribution shift. As shown in Fig. 4, both effects can result in 
overestimation of AOD. In this case, the overestimation re- 
sulted in an improvement in underestimated AOD from the 
C005 algorithm. Flowever, the standard deviation of the re- 
trieval errors is slightly lower than that of the C005 algo- 
rithm for the coarse-dominated regime, while the standard 
deviation is similar between the two algorithms for the fine- 
dominated regime (AE > 1.3). As a result, the new aerosol 
models can be inferred to mitigate systematic errors (i.e., 
MB) compared with the C005 algorithm except for severe 
dust events (AOD > 0.3, AE < 0.3). 

Figure 8 shows the scattering angle dependence of the 
retrieval errors. For the overall case of the C005 algo- 
rithm, the MB decreases gradually with increasing scat- 
tering angle and mean AOD. The decrease is a result of 
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combined effects of the systematic overestimation of low 
AOD cases (AOD < 0.2) and the underestimation of fine- 
dominated AOD with increasing AOD. In contrast, the scat- 
tering angle dependence of the MB is much lower in the 
test algorithm, indicating improved retrieval accuracy. Only 
two distinct positive peaks are present for 140° < 0 < 165°. 
For AOD > 0.3, the C005 algorithm shows systematic under- 
estimation, while the test algorithm shows positive MB for 
0 > 140° and negative MB for 0 < 140°. The underestima- 
tion in the C005 algorithm tends to worsen with increasing 
scattering angle, partially due to increasing AOD, while high 
AOD seems to correspond to low MB in the test algorithm. 

For a more detailed explanation on the scattering angle de- 
pendence of the retrieval errors, additional comparisons for 
aerosol type information are shown in Fig. 9. Note that ne- 
glecting the nonsphericity of dust particles results in underes- 
timation of AOD in the back-scattering direction, and overes- 
timation in the side-scattering direction due to the difference 
in scattering phase function between spherical and nonspher- 
ical particles. For a dust-dominated case (AE < 0.8), the data 
provided by the C005 algorithm shows an imprint of the dif- 
ference in phase functions between spherical and nonspher- 
ical particles, while the new aerosol models significantly re- 


duce the scattering angle dependence. The new aerosol mod- 
els, however, systematically overestimate the AOD regard- 
less of the scattering angle. Consequently, the small MB of 
the C005 algorithm for coarse-dominated AOD, represented 
in Fig. 7, can be explained by cancellation of the positive and 
negative errors, while the test algorithm systematically over- 
estimates AOD. For anthropogenic aerosols (AE > 1.2), nei- 
ther algorithm shows distinct features related to differences 
in the phase function, while the systematic underestimation 
of the C005 algorithm in the back-scattering direction is sig- 
nificantly reduced by the test algorithm. 

The retrieval results are further categorized with re- 
spect to the SSA in order to understand the impact of 
aerosol absorption on retrieval accuracy. Figure 10 shows 
the MB of the retrieved AODs for SSA. The AE con- 
straint is also applied to separate fine-dominated and 
coarse-dominated cases. The daily mean SSA values from 
AERONET are used in this investigation since insufficient 
data points are gathered if a time constraint is applied. 
Note that AERONET Level 2 SSA values are retrieved for 
AOD (440 nm) > 0.4. Thus, fine-dominated (AE > 1.2) and 
coarse-dominated cases (AE < 0.8) represent anthropogenic 
and dust events, respectively. For coarse-dominated case, no 
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distinct dependence on SSA is observed for either dataset, 
but systematic overestimation of dust AOD from the test 
algorithm is highlighted. However, smaller standard devia- 
tions from the test algorithm indicate systematic error, while 
C005 algorithm shows somewhat larger dispersions. On the 
other hand, relatively strong absorptivity dependence is ob- 
served in the C005 algorithm for fine-dominated case, show- 
ing increasing underestimation tendency with SSA. In fact, 
the result is confusing because no highly absorbing aerosol 
model is included in the C005 algorithm and underestimation 
of absorptivity normally results in underestimation of AOD. 
However, the complexity of error sources can produce the 
anomaly. As found in this study, both absorptivity and parti- 
cle size assumed in the retrieval algorithm can cause errors in 
the retrieved AOD. In addition, errors in the retrieved FMF 
can result in additional error in the AOD. The results from the 
test algorithm show relatively accurate AOD retrieval with 
mitigated SSA dependence. 

The sensitivity of TOA reflectance to AOD increases with 
air mass because of the increasing optical path. Thus, the 
air mass factor (AMF), defined by n?sun x m satellite where 
m — sec(((), can affect retrieval accuracy. Figure 11 shows 
the retrieval error dependence on the AMF. We expected 
the retrieval errors to decrease with the AMF because of in- 
creased sensitivity; however, the pattern is more complicated 


because the retrieval errors are functions of AOD, ©, and 
aerosol type. For both algorithms, retrieval stability increases 
(decreasing standard deviation) with the AMF partially be- 
cause of increasing sensitivity. However, the MB shows dif- 
ferent behavior; it decreases in the negative regime and then 
increases in the positive regime with increasing AMF for 
AMF < 1.6, and then gradually decreases with increasing 
AMF. For the high AOD case, however, a high AMF does not 
guarantee retrieval stability; the standard deviation is uncor- 
related with the AMF. For the test algorithm, the MB tends 
to decrease with increasing AMF except for the bifurcation 
observed for AMF < 1.5, while no dependency is observed 
for the C005 algorithm. 


6 Conclusions 

We quantitatively assessed the impact of the use of new 
aerosol models on AOD retrieval from spectral reflectance 
observed by Aqua-MODIS over the global ocean for the 
period from 2003 to 2010. AERONET inversion data and 
the optical property data of tri-axial ellipsoidal dust parti- 
cles from an existing database were used to compile AOPs 
in order to calculate LUTs, which include various aerosol 
types from absorbing to non-absorbing (0.85 < SSA < 1.00) 
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and from fine- to coarse-dominated (0.2 or 0.3 < FMF < 1.0). 
Because the MODIS C005 algorithm considers only water- 
soluble aerosols with/without humidity for fine-mode and sea 
salt/dust for coarse-mode, a noticeable difference was ob- 
served in AOD retrieval using the new algorithm with the 
consideration of various absorptivities and the size distribu- 
tion change as a function of AOD. 

Validation of the algorithms by using eight years of data 
revealed the new aerosol models to improve the AOD, with 
a regression equation of y = 0.99x +0.007 and a Pearson 
coefficient of 0.93 compared to y = 0.85x + 0.028 and 0.92 
for the C005 algorithm. The percentage of AOD data falling 
within the expected error was 64 % for the test algorithm and 
62 % for the current operational algorithm. In particular, im- 
provements were noted in the high AOD regime (AOD > 0.3) 
where the aerosol signal dominates the surface signal with a 
12 % increase in the number of reliable data points within the 
expected error. The root mean squared error (RMSE) and MB 
were also improved by the use of the new aerosol models. 

To further characterize the retrieval errors, the data were 
validated with respect to AE, scattering angle, SSA, and 
AMF. The new aerosol models mitigated the dependence 
of MB (systematic error) on the aforementioned parameters. 
However, the coarse particle-dominated AOD was still over- 
estimated and the fine particle-dominated AOD was under- 


estimated. While the systematic overestimation of the coarse 
particle-dominated AOD increased for the high AOD case, 
the results for the fine particle-dominated AOD cases were 
similar compared to the overall case. Retrieval stability, how- 
ever, was higher for the coarse particle-dominated case than 
the fine particle-dominated case, partially due to the wider 
variability in the optical properties of the fine-mode aerosols. 
In addition, the constrained analyses revealed the test algo- 
rithm to significantly reduce the scattering angle dependence 
of the retrieval error for dust-dominated cases (AE<0.8, 
AOD > 0.3), partially due to improved treatment of the non- 
sphericity of dust particles and mitigated SSA dependence 
for fine-dominated cases. The standard deviation of the re- 
trieval errors for the overall case tended to decrease with 
AMF as expected, but no distinct tendency was observed for 
AOD > 0.3. Our validation results indicate that the aerosol 
models adopted in the current MODIS operational algorithm 
need to be updated to achieve better accuracy. Further anal- 
yses on FMF and SSA retrievals are required to obtain a 
better understanding of the various error sources contribut- 
ing to AOD retrieval and corresponding improvements in the 
aerosol retrieval algorithms. 
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